Modeling the actinides with disordered local moments 
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A first-principles disordered local moment (DLM) picture 
within the local-spin-density and coherent potential approx- 
imations (LSDA+CPA) of the actinides is presented. The 
parameter free theory gives an accurate description of bond 
lengths and bulk modulus. The case of S-Pu is studied in 
particular and the calculated density of states is compared to 
data from photo-electron spectroscopy. The relation between 
the DLM description, the dynamical mean field approach and 
spin-polarized magnetically ordered modeling is discussed. 



I. INTRODUCTION 

The elemental actinide metals, Pu in particular, ex- 
hibit several unique features. They are among the most 
complex elements in nature, with a rich set of allotropes 
[1], of which several have complex low symmetry crys- 
tal structures [2]. Pu, for example, is the only element 
with seven condensed matter phases at zero pressure of 
which one (the <5-phase) demonstrates negative thermal 
expansion. The understanding of these anomalous prop- 
erties is a serious challenge. Several of the main features 
of the actinides can be understood in terms of the pro- 
gressive filling of the 5f shell. In the light actinides, Th 
to Np, the 5f electrons are itinerant and participate in 
the bonding, whereas in the heavier actinides, Am to Cf, 
the 5f electrons are localized and exhibit behavior more 
similar to the lanthanides. Plutonium takes a particu- 
lar place between these extremes. The low temperature 
cv-phase has been shown to be well described with itiner- 
ant bonding 5f electrons whereas the high temperatures 
phases with their increased volumes suggests a localized 
or partly localized configuration. This behavior indicates 
a Mott-like transition of the f-electrons similar to the a-j 
transition in cerium [3]. 

There have been several first principles approaches to 
calculating the properties of the actinides [2,4-8,10-15]. 
Although the light actinides are well described within 
density functional theory, the local density approxima- 
tion (LDA) fails in general to describe the 5f localiza- 
tion which occurs in actinide compounds, elemental Pu at 
high temperature, and in the actinides past Pu. A num- 
ber of techniques have therefore been developed to in- 



clude correlation effects beyond LDA in order to describe 
the partially localized 5f electronic structure. These in- 
clude LDA+U [16,17], orbital polarization [15], and the 
Mixed Level Model [13]. Of particular relevance to the 
work presented here are recent attempts to combine LDA 
and Dynamical Mean-Field Theory (DMFT). DMFT has 
been shown to be a very useful approximation at the 
consideration of strongly correlated electron systems (for 
review, see [18]). LDA combined with the DMFT may 
provide a first principles technique to study correlated 
electron materials, and there have been several attempts 
to apply implementations of LDA+DMFT [19,20]; for 
a review, see [22,23]. Recently this LDA+DMFT ap- 
proach has been applied to the Pu problem [21]. The 
LDA+DMFT approach gives an opportunity to describe 
correlation effects on the electronic structure and proper- 
ties of d- and /- electron systems. However, the technique 
is cumbersome and it is not completely ab initio because 
of the problem with the choice of U (see, e.g. [22]). In 
addition, implementations of DMFT are computation- 
ally intensive which makes the calculation of complex 
structures difficult. Sometimes, because of the complex- 
ity of the calculations, uncontrollable approximations are 
made, e.g., using a single U parameter instead of a com- 
plete interaction matrix [24,25,21]. 

It is commonly accepted that accounting for Hubbard 
correlation effects is of crucial importance for /-electron 
systems. On the other hand, some of these effects can be 
taken into account, in an approximate way, in the frame- 
work of more traditional density functional techniques. 
This is the approach taken in the present study where we 
have modeled the actinides by means of Disordered Local 
Moments (DLM) [26] within the Local Spin Density Ap- 
proximation and the Coherent Potential Approximation 
[27] (LSDA+CPA). The purpose of the present work is 
to investigate and possibly demonstrate that correlation 
effects beyond the standard local density approximation 
can be simulated, at least partially, by means of a pa- 
rameter free first principles DLM approach, and to see 
if this picture gives an adequate description of the ac- 
tinides. The DLM picture, even if it is insufficient for a 
complete description, might lead to some insights for the 
understanding of the electronic structure of the actinides. 

Consider, first, the status of the DLM approach in the 
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many-body lattice models like Hubbard or s-f exchange 
("Kondo lattice" ) model. The key point is an equivalence 
between a many-body interacting system with Coulomb 
on-site interactions and a one-electron system in fluctu- 
ating charge and spin fields. This equivalence, which can 
be proven by the Hubbard-Stratonovich transformation, 
is a base of spin-fluctuation theories of itinerant-electron 
magnetism [28]. In a complete theory, the charge and 
spin fields are dynamically fluctuating both in space and 
time. However, a "static approximation" , where we ne- 
glect the dynamics of the fluctuations, [29,30] captures 
an important part of the correlations while greatly sim- 
plifying the formalism, and may be sufficient for many 
problems of interest. In this case the correlated sys- 
tem is described in terms of a DLM alloy and a CPA 
for this alloy becomes equivalent to the "Hubbard III" 
approximation [31] for the original many-body problem 
(see [32]). The Hubbard III approximation, and there- 
fore also the corresponding DLM description, gives qual- 
itatively the correct picture of the electronic structure 
both in atomic and broad-band limits and can describe 
the metal-insulator (Mott) transition for half-filled bands 
when U is of the order of the bandwidth [31]. The most 
important difference with the more elaborate DMFT pic- 
ture [18] is an absence of the "Kondo peak" near the 
Fermi level at the metallic side of the transition. How- 
ever, this Kondo peak has a small spectral weight and, 
generally speaking, it is not obvious that it is of cru- 
cial importance for the description of total energy and 
related characteristics. It is shown in Ref. [33] that the 
Hubbard III approximation can be rigorously justified, 
in the s-f exchange model, in the limits of large space 
dimensionality and classical spins (S — > oo) whereas for 
DMFT only the first limit is essential. Note, that the 
Kondo resonance is essential also for a correct description 
of the properties of electron states in doped Mott insu- 
lators (large -U limit near the half filling). On the other 
hand, for half filling the Hartrcc-Fock-likc descriptions of 
the band-split states arc, in general, quite adequate [34]. 
Physically, this is the consequence of a freedom in the de- 
scription of occupied bands in terms of Bloch functions, as 
well as Wannier functions or any other orthogonal linear 
combinations. The spatial fluctuations of the exchange 
on-site field can lead to a splitting of the energy spec- 
tra (provided the fluctuations are larger than the band- 
width) which corresponds, in the Hubbard model, to the 
Hubbard band splitting. The self-energy in the DLM 
picture is energy dependent and has an imaginary part 
describing the damping of the electrons on spin fluctu- 
ations. This distinguishes the DLM approach not only 
from magnetically ordered LSDA calculations but also 
from LDA+U, SIC, and the Hartree-Fock approximation. 
In this sense, the DLM-CPA approach can be considered 
as a particular case of the "LDA++" approach [20] with 
a local, energy dependent, complex self-energy. However, 
in contrast to schemes taking into account Hubbard cor- 



relations, it is completely ab initio. A shortcoming of 
this approximation is an incorrect description of electron 
damping near the Fermi level; i.e., the absence of the 
Kondo resonance, which is a consequence of the quantum 
character of spin, as well as problems describing localiza- 
tion in systems that do not have half filled (or completely 
filled) electron shells. However, this shortcoming may be 
of minor importance for the description of electron en- 
ergy spectrum at large energy scales as well as for the 
calculations of total energy and related characteristics 
such as the equilibrium volume and elastic moduli. Note 
also that this finite damping, because of electron scatter- 
ing by spin fluctuations, is a physically correct picture 
for high enough temperatures. 

Comparing the many-body lattice models with the 
density functional approach one can face with the well- 
known "Hubbard U vs. Stoner J" problem, i.e., with the 
inadequacy of the LDA approach near the atomic limit 
[16,20]. On the other hand, for moderately correlated 
systems such as, e.g., 3d metals, the main correlation 
effects are connected with spin degrees of freedom and 
can be described, in principle, basing on the LSDA elec- 
tronic structure. Of course, it is difficult to say a priori 
where the boundary is between 'moderately correlated" 
and "strongly correlated" systems. We will demonstrate 
that, at least for early actinides a LSDA-DLM descrip- 
tion of the electronic structure turns out to be rather 
successful. 

Due to the localized character of the spin moments, 
any ordered magnetic structure will resemble the re- 
sult of the DLM description. It has often been noted 
that localization effects could partly be accounted for by 
way of exchange spin splitting of the 5f band, and de- 
spite an unphysical long range magnetic ordering (absent 
in the DLM approach), several spin-polarized ferromag- 
netic (FM) and anti-ferromagnetic (AFM) calculations 
have been performed [4,6,10,14,15]. The presented DLM 
model gives a natural generalization to the paramagnetic 
state. 

First we discuss some calculational details and then 
present the calculated volumes and bulk moduli for the 
actinide series. We then look at the results for <5-Pu in 
detail, and the results are compared with different alter- 
native approaches and experiments. 

II. CALCULATIONAL DETAILS 

All total energies and densities have been calculated 
self-consistently within the framework of density func- 
tional theory [35,36], in the local-density approxima- 
tion (LDA) in the non-magnetic cases, and within the 
local-spin-density approximation (LSDA) for the spin- 
polarized systems [37-39], with the local exchange- 
correlation functional by Perdew et al. [40]. 
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We used the basis set of the s, p, d, and / linear 
muffin-tin orbitals (LMTO) in the tight- binding repre- 
sentation and the atomic sphere approximation (ASA) 
for the crystal potential [41-44]. The method was im- 
plemented within the scalar relativistic Green's function 
technique [45-48]. Spin-orbit coupling was not included. 
The disorder of local spin moments was treated within 
the coherent potential approximation (CPA), and other 
details relevant for the present calculation can be found 
in Rcfs. [47,49]. 

Since the interatomic distance is known to be the most 
important factor that determines the energetic of an ac- 
tinide and since we are mainly interested in trends rather 
than in a detailed quantitative description of all the 
different phases of the actinides, calculations were per- 
formed for the fee crystal structure only. 

III. DISORDERED LOCAL SPIN MOMENT 
PICTURE OF THE ACTINIDES 

Figure 1 shows a comparison between experimental 
and calculated Wigner-Seitz radii for the actinide met- 
als. For the non-magnetic LDA calculation, we find a 
parabolic behavior typical of bonding through a series 
such as the transition metals. This is in full agree- 
ment with the experimental values for the early actinides, 
but fails do describe the volumes of the heavy actinides. 
Within the DLM picture on the other hand, we find a 
substantial improvement. For the light actinides, Th-U, 
results are identical to the non-magnetic LDA calcula- 
tion, though for Np we find a slightly increased volume. 
For the later actinides we find an abrupt volume increase 
in close agreement with the experimental values with the 
transition taking place at Pu, which has an intermediate 
volume close to the value of the high temperature <5-phase 
in the fee structure. The results in Fig.l are quite close 
to the data presented by Skriver et al. [4]. 

We find similar results for the bulk moduli. In fig- 
ure 2 we display a comparison between the experimental 
bulk modulus for the ground state structures and the 
calculated bulk modulus in the fee structure. The LDA 
calculations gives a fairly good description of the early 
actinides but fails completely for the later actinides. The 
DLM picture on the other hand gives a significantly im- 
proved result where, for example, the bulk modulus of 
Pu is reduced by a factor of two, though still somewhat 
higher than the experimental result for 5-Pu. 

Given the fact that the results are from parameter free 
ab initio calculations without considering the exact crys- 
tal structures and spin-orbit coupling the agreement be- 
tween the calculated equilibrium volumes and bulk mod- 
uli with experimental values is remarkably good and the 
results clearly indicate the ability of the LSDA+CPA 
approach within the DLM picture to model important 



correlation effects beyond LDA in the 5f band, without 
incorporating long ranged magnetic ordering. 

The mechanism behind the improved description of the 
equilibrium volume and bulk modulus is the formation 
of local disordered moments. This gives rise to a band 
splitting and partial localization of the 5f electrons which 
reduces the bonding resulting in an increased lattice con- 
stant and reduced bulk modulus. For Th and Pa local 
spin moments are still quenched, and for U only a negligi- 
ble moment of about 0.07 /xb is formed. Within our DLM 
scheme, neptunium has a moment of 1.69 /zb- Though 
this has little effect on the equilibrium lattice constant, 
the bulk modulus is reduced by a factor of two. For 
Pu a spin moment of 4.76 /j,b is found in the DLM cal- 
culation, which corresponds to an almost complete spin 
polarization of the 5f electrons. For Am and Cm a mo- 
ment of 6.57 and 6.90 /is is formed, respectively. For Bk 
the disordered local moment is slightly reduced to 5.61 
/ie- Notice that fully relativistic calculations, including 
the spin-orbit interaction, would reduce the moments. In 
the case of 5-Pu from about 5 /is to less than 2 ^b [6] 
and for Am the total moment can be expected to be ~ 5 
/xb [7]. However, the formation of a J=0 atomic ground 
state comes natural from a completely localized 5f shell, 
[5] and is consistent with experimental data. 

IV. 5-PU 

Figure 3 shows the calculated density of states (DOS) 
for Pu in comparison with the results of photo-electron 
spectroscopy [50]. The LDA curve has qualitatively a 
very different behavior compared to experiments. For ex- 
ample, the peak at the Fermi level seen in the photoemis- 
sion spectra is absent. The DLM DOS on the other hand, 
is in qualitative agreement with experimental data. The 
DOS is similar to the LDA-DMFT spectra by Savrasov 
et al. [21] (calculated at 600 K), but the quantitative 
agreement for the DLM picture is considerably better. 

The electronic structure is strongly modified in the 
DLM picture compared to the LDA calculation. The 5f 
band is split, as shown in Fig. 4, due to the formation of 
local disordered moments, leading to an effective energy 
lowering of the occupied 5f band mass. The same effect is 
found also for the AFM calculation. This leads to a de- 
creased bonding, increased equilibrium lattice constant 
and a reduced bulk modulus. The effect is similar to the 
results of LDA+U, LDA++, and LDA-DMFT. In con- 
trast to some of these models we here have an improved 
agreement with the photo-electron spectra and the result 
is similar to what was found in the Mixed Level Model 
by Eriksson et al. [13]. However, though a direct full 
comparison between the Kohn-Sham eigenspectra and 
the photo-emission data may not be possible, the data 
clearly indicates a resonance at the Fermi level which is 
not present in the LDA description. 
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In Fig. 3 we also display the density of states for an an- 
tiferromagnctic (AFM-LSDA) calculation. Though, re- 
sembling the general DLM 5f band split, the DOS is in 
less agreement with photoemission data. 

V. DISCUSSIONS 

There are two principle approaches for describing the 
a — 7 transition in Ce and the a — 6 transition in Pu: the 
Mott transition model by Johansson [3] (for Pu, it has 
been considered in Refs. [9,10]) and the Kondo-collapsc 
model of Allen and Martin [51]. In the context of DMFT, 
the appearance of the Kondo resonance is a feature of 
the Mott transition. An important question is thus what 
is more important for the energetics of the transitions: 
the appearance of the Hubbard gap or the formation of 
the Kondo resonance? The first feature is taken into 
account in the DLM scheme whereas the second one is 
not. It was already mentioned that the Kondo resonance 
has small spectral weight, so it is not obvious that ig- 
noring it is an essential shortcoming. Moreover, both 
for Ce and Pu there are some arguments connecting the 
Mott transition to the peculiarities of the atomic elec- 
tronic structure (atomic collapse [52]). These peculiar- 
ities can obviously be taking into account in the local 
spin-density approach and thus in the DLM description. 
Note that this collapse phenomenon, leading to a sharp 
dependence on the band width and lattice constant, leads 
to a first-order phase transition. Therefore the differ- 
ence between the description of the Mott transition in 
the DLM or the Hubbard III approximation (continuous 
transition [29,33]) and DMFT (first-order transition, see 
Ref. [53] and Refs. therein) which is very essential for the 
Hubbard model, does not play a serious role for real ac- 
tinide systems where the transition is essentially of first 
order for "quasi-one-electron" reasons. 

As for the description of magnetic properties at high 
temperatures T, the static approximation gives a quali- 
tatively correct picture. The magnetic susceptibility \ is 
proportional to < e 2 > /T, where e is the exchange on- 
site field. The susceptibility corresponds to the paramag- 
netic Pauli spin susceptibility provided that the exchange 
splitting fluctuations are Gaussian near the point e = 0, 
and to the Curie- Weiss susceptibility provided that spon- 
taneous spin splitting exists [54]. The description of finite 
temperature magnetism of transition metals in the DLM- 
CPA approach [26] is different from the DMFT treatment 
[55] only by the consideration of spins in a classical way; 
in the Heisenberg model, it corresponds to the appear- 
ance of the multiplier S 2 instead of the correct factor 
S(S + 1) in the Curie constant. 

The DLM picture leads to an expected Curie- Weiss 
behavior of the magnetic susceptibility for 5-Pu. How- 
ever, recent experiments with Al and Ga stabilized S-Pu 
shows only a vannishingly small temperature dependence 



of the magnetic susceptibility, see Fig. 2 in Ref. [57] . This 
is also consistent with other susceptibility measurements 
[58,59]. The possible reason, explaining the absent Curie- 
Weiss behavior, could cither be due to crystal field split- 
ting and the formation of a non-magnetic multiplet [13] 
or because of the occurrence of a Kondo resonance that 
may diminish the temperature dependence of the suscep- 
tibility provided that the resonance width ( "Kondo tem- 
perature" ) is larger than T. At least the first explanation 
could possibly be included separately into the DLM pic- 
ture, but the second one can be described only within the 
framework of DMFT. At the same time, it is important 
to realize that the intention with the DLM picture is to 
model some of the main characteristics of the energetics 
of the actinides, and it does not necessarily describe the 
magnetic properties correctly. 

VI. SUMMARY 

In summary we have presented a first principles Dis- 
ordered Local Moment (DLM) method within the local 
density approximation with the disorder treated within 
the coherent potential approximation. The DLM picture 
gives an reasonably good description of bond lengths and 
bulk modulus for the actinide series. The equivalence 
between the DLM picture and the Hubbard III approxi- 
mation and their relation to the DMFT description was 
discussed and it was argued for that the DLM picture is 
related to DMFT through a static approximation. The 
DLM density of states compares well with photoemission 
on <5-Pu, in contrast to that obtained from the LDA or 
the magnetically ordered AFM configuration. In general, 
it is found that the DLM picture gives a considerable im- 
provement over the LDA results and quantitatively good 
agreement with experiments. 
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FIG. 1. 

Comparison between the experimental [5] and calcu- 
lated atomic Wigner-Seitz radius, R\vs> for the actinide 
metals. 
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FIG. 2. 

Comparison between the experimental and calculated 
bulk modulus for the actinide metals. The experimental 
data are given in a) Ref. [5] and b) [56]. 

FIG. 3. 

Comparison between the calculated DOS within the 
DLM description, non-magnetic LDA, an antifcrromag- 
netic (AFM) structure, and the photo-emission spectra 
(PES) [50] for 5-Pu. 

FIG. 4. 

Calculated DOS for <5-Pu within the DLM description 
in comparison with the non-magnetic LDA DOS and an 
antiferromagnetic (AFM) LSDA DOS. 
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